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THEORY  OF  ELECTRON  CYCLOTRON  MASER 
INTERACTIONS  IN  A CAVITY  AT  THE  FUNDAMENTAL 
AND  HIGHER  CYCLOTRON  HARMONICS 

I 

L INTRODUCTION 


I i 


V ! 

I 


In  recent  years,  there  has  been  increasing  interest  in  an  electromagnetic  radiation 
mechanism  known  as  the  electron  cyclotron  maser.1  ~3  In  addition  to  its  fundamental  impor- 
tance as  a new  scientific  phenomenon,  it  has  been  the  basis  for  a new  type  of  microwave 
devices4,5  (gyrotrons)  capable  of  generating  microwave  at  unprecedented  power  levels  in  the 
centimeter  through  submillimeter  wavelength  regime.  This  new  source  of  high  power  mi- 
crowaves has  shown  gteat  promises  for  plasma  heating6  in  fusion  devices  and  opened  new 
areas  of  radar  applications.  The  cyclotron  maser  radiation  mechanism  originates  from  a rela- 
tivistic effect.  For  the  simplicity  of  illustration,  we  consider  an  ensemble  of  monoenergetic 
electrons  in  the  presence  of  a uniform  magnetic  field  B0  and  assume  that  the  electrons  have 
no  velcoity  along  the  magnetic  field,  hence  they  all  gyrate  with  the  same  Larmor  radius.  Initial- 
ly, the  phases  of  electrons  in  their  cyclotron  orbits  are  random,  so  no  radiation  will  be  emitted. 
But  phase  bunching  can  occur  because  of  the  dependence  of  the  electron  cyclotron  frequency 
on  the  relativistic  electron  mass.  Those  electrons  that  are  decelerated  in  the  wave  electric  field 
become  lighter,  rotate  faster,  and  hence  accumulate  phase  lead  while  those  electrons  that  are 
accelerated  rotate  slower  and  accumulate  phase  lag.  This  will  result  in  phase  bunching  such 
that  the  electrons  radiate  coherently  at  the  frequency 

<o=rs{le/y 
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where  - eB^fmc,  y is  the  relativistic  factor,  and  s is  an  integer.  Original  investigations  of 
the  physical  processes  were  carried  out  in  the  late  1950’s  by  Twiss, 1 Schneider,2  and 
Gaponov,3  and  early  experimental  observations  were  made  by  Hirshfield  and  Wachtel,7  and 
Chou  and  Pantell.8  Further  theoretical 9 ~28  and  experimental 29  ~42  studies  as  well  as  review 
works4,5  have  appeared  frequently  in  literature.  The  azimuthal  bunching  mechanism  described 
here  and  the  well  known  Lorentz  force  induced  axial  bunching  mechanism  are  evidently  of 
different  physical  origins;  however,  as  shown  in  a recent  comparative  study,15  the  two  mechan- 


isms always  combine  in  such  a way  as  to  offset  one  another.  The  azimuthal  bunching  mechan- 
ism dominates  for  waves  with  phase  velocity  greater  than  the  speed  of  light  and  vice  versa, 
hence  the  cyclotron  maser  radiation  can  be  generated  only  in  a fast  wave  structure. 


Theoretical  cyclotron  maser  studies  have  generally  been  based  on  two  models,  each 
corresponding  to  a common  experimental  configuration.  In  the  first  model,9-16  the  elec- 
tromagentic  wave  grows  as  an  instability  due  to  the  velocity  space  anisotropy  of  the  electron 
medium.  It  applies  to  travelling  wave  amplification  in  waveguide  structures.  In  the  second 
model,17  -28  on  which  the  present  analysis  will  be  based,  the  electron  beam  interacts  with  the 
constant  amplitude  standing  wave  of  a cavity  structure.  It  applies  to  beam  sustained  oscillations 
in  a finite  Q cavity. 

In  previous  analytical  studies  17,18,21  -24,27  of  the  second  model,  it  has  been  generally  as- 
sumed that  the  electron  Larmor  radius  is  small  compared  with  the  radial  scale  length  of  the 
field  structure  so  that  field  variations  in  the  radial  direction  can  be  neglected.  This  assumption 
simplifies  the  analysis  considerably  while  still  leading  to  accurate  results  for  the  fundamental 
cyclotron  harmonic  interaction.  However,  information  concerning  higher  cyclotron  harmonic 
interactions  has  been  precluded  under  this  assumption  because  the  higher  harmonic  interac- 
tion is  essentially  a finite  Larmor  radius  effect.  A unique  advantage  of  the  cyclotron  maser 

comes  from  its  capability  to  generate  high  power  submillimeter  waves.  Using  the  relation 
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o>  = sile/y  [Eq.  (1)],  we  find  that  generation  of  submillimeter  waves  through  the  fundamental 
cyclotron  harmonic  (s  — 1)  requires  an  external  magnetic  field  in  excess  of  100  kG,  apparently 
too  high  to  be  practical.  Thus  interactions  at  high  cyclotron  harmonics  become  especially  ad- 
vantageous if  large  magnetic  field  is  to  be  avoided.  In  view  of  these  considerations,  our  em- 
phasis in  the  present  study  will  be  on  the  higher  cyclotron  harmonic  interactions. 

Starting  from  the  relativistic  Vlasov  equation  and  the  Maxwell  equations,  we  calculate  the 
linear  responses  of  an  annular  electron  beam  as  it  passes  through  the  constant  amplitude  fields 
of  a cavity  structure.  In  contrast  to  previous  treatments,  the  exact  spatial  field  variations  (radial 
variation  in  particular)  have  been  incorporated  and  the  electron  Larmor  radius  has  been  kept 
arbitary.  From  the  linear  beam  responses,  the  beam-wave  coupling  coefficient  and  the  beam 
energy  gain  function  are  derived  for  harmonic  interactions.  On  the  basis  of  these  results,  con- 
ditions for  maximum  beam-wave  coupling  with  respect  to  such  parameter  as  beam  position, 
beam  energy,  and  cavity  length  can  be  found.  It  is  shown  that,  for  each  mode  of  operation, 
there  exists  a threshold  beam  power  below  which  the  cavity  oscillations  can  not  be  started.  In 
general,  the  higher  the  cyclotron  harmonic,  the  higher  the  thresold  beam  power,  and  the 
longer  the  cavity,  the  lower  the  threshold  beam  power.  It  is  found  that  such  threshold  condi- 
tions can  become  very  restrictive  for  higher  cyclotron  harmonics,  hence  in  such  operations  it  is 
desirable  to  have  parameter  fully  optimized  so  as  to  lower  the  threshold  conditions.  In  this  re- 
gard, the  analyses  to  be  presented  could  be  of  considerable  importance  in  assessing  the  feasibil- 
ity of  experimental  goals  and  in  selecting  the  proper  modes  and  parameters  to  achieve  them. 

In  Section  I,  the  problem  is  formulated  for  a general  beam  distribution  function.  In  Sec- 
tion II,  we  specialize  to  a distribution  function  consistent  with  the  beam  generated  in  a 
magnetron-type  electron  gun  commonly  used  for  cyclotron  maser  experiments.  In  Section  III, 
the  final  results  are  presented,  together  with  physical  interpretations  and  numerical  examples. 

A comparison  is  made  between  the  nature  of  cyclotron  maser  interactions  in  a cavity  structure 
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and  that  in  a waveguide  structure.  A numerical  code  has  been  developed  as  an  independent 
check  of  the  analytical  results.  Section  IV  discusses  the  validity  of  the  approximations  used  and 
some  nonlinear  aspects  of  the  problem. 
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n.  MODEL  AND  FORMULATION 

Figures  la  and  b show  the  configuration  of  the  electron  cyclotron  maser  system  under 
study.  It  consists  of  an  annular  electron  beam  propagating  inside  a circular  cross-section  cavity 
(radius  rw,  length  L).  The  axis  of  the  electron  beam  coincides  with  that  of  the  cavity.  The 
electrons,  guided  by  an  applied  uniform  magnetic  field  (£oe2).  move  along  helical  trajectories. 
The  electrons  have  a substantial  part  of  their  kinetic  energy  in  the  form  of  transverse  gyromo- 
tion  and  the  rest  in  the  form  of  axial  motion.  Two  types  of  electron  orbits  are  implicit  in  this 
model.  The  more  common  type  is  shown  in  Fig.  la  where  the  electron  orbit  does  not  encircle 
the  axis.  The  other  type  has  the  electron  orbit  encircling  the  axis. 

The  following  simplifying  assumptions  are  made: 

1.  The  electron  beam  distribution  function  and  the  cavity  fields  are  both  independent  of 
the  azimuthal  angle  6. 

2.  The  beam  is  sufficiently  tenuous  so  that  its  self  electrostatic  and  magnetic  fields  are 
negligible  compared  with  the  cavity  fields. 

3.  The  cavity  fields  are  of  first  order  compared  with  the  applied  magnetic  field  B0\  and 

4.  The  perturbed  distribution  function  /(1)  is  of  first  order  compared  with  the  initial  dis- 
tribution function  /0. 

For  cyclotron  maser  interactions,  it  is  well  known 2 3  4 that  the  beam  couples  much  more 
strongly  with  the  TE  mode  than  the  TM  mode.  This  together  with  assumption  1 will  thus  res- 
trict our  consideration  to  the  TE0nl  cavity  mode,  where  n and  / are  the  radial  and  axial  eigen- 
mode numbers,  respectively.  The  field  components  of  the  TE0nj  cavity  mode  are 
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and 


£e(1  3 — Eqq  Jx  ( kn  r)  sin  kzz  cos  «f, 
£r(1)  — (kzc/u>)  EgQ  Jx  ( k„r ) cos  kzz  sin  tut, 


(2) 

(3) 


B/1 * — —(knc/o> ) £90  70  (x„r)  sin  fcz2  sin  mt, 

(4) 

where  kz  = ttI/L,  kn  = x,,//^  x„  is  the  n-th  nonvanishing  root  of  Jx  (x)  — 0 , 
on  = (kz  + /c2)1/2c  is  the  wave  frequency,  and  superscript  "(1)"  indicates  first  order  quanti- 
ties. 


1 


The  dynamics  of  the  electron  beam  are  described  by  the  realtivistic  Vlasov  equation, 

+ v • — e(E  + — v x B)  • =0. 

dr  dr  c dp 


(5) 


Linearizing  Eq.  (5)  by  the  use  of  the  ordering  schemes  of  assumptions  3 and  4,  and  in- 
tegrating the  resulting  equation  by  the  method  of  characteristics,  we  obtain 


/(1)  (p,  r,  t)  - f dr'e[E(1)(r',  r')  +-v' x B(1)(r ',/')]  ~/0(r’,p'), 

'-*/»«  c dp  (6) 

where  v — p/ym,  y = [1  + (/r2  + p2  )/m2c2] 1/2,  the  primed  quantities  are  treated  as  func- 
tions of  1',  and  r,  p are,  respectively,  the  values  of  r',  p'  at  t‘  — t.  The  integration  over  t'  is  to 
be  carried  out  along  the  unperturbed  (helical)  orbits  of  the  electrons.  Along  its  orbit,  an  elec- 
tron feels  the  axial  as  well  as  the  radial  variations  of  the  electromagnetic  fields,  E (1  * and  B (1  \ 
given  by  Eqs.  (2)  through  (4).  An  electron  located  at  the  axial  position  z at  time  t enters  the 
cavity  at  time  t —z/vz.  This  determines  the  limits  of  the  t -integration  in  Eq.  (6). 


Methods  for  evaluating  the  integral  in  Eq.  (6)  are  standard43  and,  for  the  present  prob- 
lem, involve  the  use  of  the  following  Bessel  function  identity44: 


e 


is0t 


Js  (X,  ) 


“ i Js+,'  (x2)  Js  (x3)  e"'*2, 

j'  — -00 


(7) 
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where  xl,x2,  x3,  0j  and  02  are  related  through  the  triangle  shown  in  Fig.  2 The  lengthy  alge- 
bra will  not  be  detailed  here.  Instead,  in  Sec.  IV,  we  will  present  detailed  physical  interpreta- 
tions and  an  independent  numerical  check  of  the  final  results. 

The  perturbed  distribution  function  /(1)  as  solved  from  Eq.  (6)  assumes  the  following 


f(  l)  _y(l)  + f 0) 


where 


L. 


fd)  — eikiz  - «•"  _ Pz  | Vo  + px.  Vq  1 

+ I 2(u  kzc  ymc  J bpL_  ymc  dpz  J 


I I 

5 — — oo  j'  a — oo 


Js'(k„r)  Gss.(knrL)Xtxp[-is'(<t>  -0)] 

I S I ■ — 

oi  ~kz\z  — sil  e/y 


and  /IP  is  given  by  Eq.  (9)  with  w replaced  by  — w.  In  Eq.  (9),  /m{Z)  indicates  the  imaginary 
part  of  Z,  Gjj-Or)  = Js+s(x)dJs  ( x)/dx , flf  = eB^/mc,  rL  = pjmile  <f>  and  0 are  respective- 
ly, the  polar  angles  of  the  momentum  and  position  vectors  ( see  Fig.  3),  and 
X = 1 — exp  [/(to  ~kz\z  — s(le/y  )z/\z).  The  perturbed  azimuthal  current,  7e(1),  can  be 
written  in  terms  of /*’  * as 


oo  2 n 


jtt)  _ -e  J f^x)  \ed3p  - —e  f p±dp±  f dpz  j d4> /(1 ) v±  sin  0, 


where  <t>  = <t>  — 0 . 


Inserting  Eq.  (8)  into  Eq.  (10)  and  carrying  out  integrations  by  parts  over  pk  and  pz , we 


obtain 


where 


i,u> 
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0.1  r - . oo  oo 

£ I JAk„r) 

4 fTlCj)  j a —co  5'  ■ — 00 


Jpa  j U/o  4 --^,+n‘l 

0 * —00  0 

•|gsj.  (knrL)p*_  l(w2  - k*c2)X/ym2c2  + /z«  (1  - JO 


• (w2  - k2c2  - kzc2c/\z)  /mc2pz]  y 2«  2 


-1  < -1 


- (w  - kzpzlym)X[2G„(knrL)  + Vi  Gss  (knrL )]  y * 


e = oj  - kzpzlym  - sile/y. 

Gjj  (x)  - dGss  (x)/dx. 

and  jfl2  is  given  by  Eq.  (12)  with  &>  replaced  by  -o. 

The  time  averaged  power  gain  (P)  of  all  the  electrons  in  the  cavity  is  then 

2tr 

<0  rw  L 


P - a>  jf  dt  f rdr  f dz  Jgix  ( £#°  * . 


Equation  (13)  is  a general  expression  for  the  beam  power  gain.  In  the  following  section, 
we  specialize  to  a particular  distribution  function  which  is  considered  to  be  most  ideal  for  cy- 


clotron maser  operations. 


III.  DISTRIBUTION  FUNCTION  FOR  A MAGNETRON  TYPE  ELECTRON  BEAM 


Before  Eq.  (13)  can  be  evaluated,  a proper  distribution  function  for  the  electron  beam  has 

to  be  constructed  from  the  constants  of  motion  of  the  system:  P±  PZ,  and  Pg,  where 


a \ £ * 

Pg  - ym  rpL  sin  4>  - — — Aq  r . 
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We  are  interested  in  a distribution  function  in  which  the  electrons  have  arbitrary  spread 
in  pL  and  pz  , with  their  guiding  centers  uniformly  distributed  on  a cyclinder  of  radius  r0  (see 
Fig.  la).  Such  a distribution  function  is  represented  by 

/0  - K 8 (rl  - 2cPe/eB0  - r§)  g {px,pz), 

(15) 

where  8(x)  is  the  Dirac  delta  function,  g(px,pz)  is  an  arbitrary  function  of  p±  and  p,  (to  be 
specified  later),  and  K is  a normalization  constant.  Note  that  in  Eq.  (15)  rL  « pjmtl  e is  also 
a constant  of  motion,  and  K is  to  be  determined  from  the  condition, 

f fa  2irrdrdip  — N, 

J (16) 
where  N is  the  line  density  of  the  electron  beam  (i.e.  number  of  electrons  per  unit  length). 


After  some  algebra,  Eq.  (15)  can  be  written 


/0  — Mr  _1  (8  (<6  — ) + 6 (<f>  ~ rr  + <£o  ) 1 £ ~ rO  S (r^ 

x [ (r2  -rf)  (r|  -r2)]  _1/2 


where 


<60  -sin  1 ((r2  + rl  - r$)/2rrL], 


1,  for  x > 0, 
0,  for  x < 0, 


r\  “ ko  ~ ri}‘ 


r2  " r0  + rL- 


and  g (px,pz)  satisfies  the  normalization  condition. 


oo  oo 

2nf  px  dpx f dpzg  (jjvpz)  - 1. 

0 * —oo 


Equation  (17)  specifies  an  electron  beam  in  which  all  electrons  have  the  same  guiding 

center  position  (r0).  Through  a simple  super  position  procedure  (e.g.  specifying  a weighing 
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function  of  r0  and  carrying  out  an  integration  over  r0),  Eq.  (17)  can  be  readily  extended  to 
treat  a beam  with  a spread  in  electron  guiding  center  positions.  In  the  remaining  part  of  this 
paper,  however,  we  shall  concentrate  on  a cold  monoenergetic  beam  with  no  spread  in  guiding 
center  positions.  Thus,  we  let 


g (PX.P2)  =8  (pa..-pifl)8  (pz  - pzo)/2npx. 


(18) 


where  p±0  and  pz0  are,  respectively,  the  perpendicular  and  parallel  momenta  of  the  electrons. 


Equation  (17)  together  with  Eq.  (18)  is  an  idealized  representation  of  the  magnetron-type 
electron  beam  commonly  used  in  cyclotron  maser  experiments35  ~37,41.  From  the  analytical 
point  of  view,  such  a distribution  function  has  some  advantages.  First,  it  leads  to  results  which 
can  be  physically  interpreted.  Second,  it  gives  simple  analytical  formulas  for  the  determination 
of  optimum  beam  positions  and  energies.  Finally,  it  serves  as  the  oasis  of  super  position  for 
forming  a beam  with  arbitrary  thermal  and  guiding  center  spreads.  One  notes  that  the  method 
of  super  position  is  applicable  here  because  the  self  fields  of  the  beam  have  been  neglected. 
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IV.  RESULT 


A.  Beam  power  gain 

Combining  Eqs.  (2),  (11),  (12),  (13),  (17)  and  (18),  we  obtain  an  explicit  expression  for 
the  time  averaged  beam  power  gain. 


P - Px  + P2  + /»3  + P4 


where 


(19) 


P\  - vc2E^T(Sy0(o)  1 £ ( ~ ffs(k„r0,  k„rL)r0*oA  3 { (to2  - k2c2)L 

• [2  (1  — cos  A ) — A sin  A ] + kzc2 A (A  sin  A + cos  A — 1 )) 

+ Os(Vo>  knrL>L*  _2(<u  ~*;v« 0)(1  "COS  A)), 

(20) 

P'1  - *-t:2££)T(8y0to)  _l  £ (#s  (Vo- Vz-M-ioA  _2A  _1{  [:(a> 2 -k2c2)L 

5 — - oo 

• ( (1  — cos  A ’)  (1  + AA  _1 ) — A sin  A ] + kzc2 A [A  sin  A — AA  _l  (1  — cos  A ’)  ]} 

- Qs(k„r0,  k„rL)L\  -1A  _1  (to  -/czvz0)(l  - cos  A )]. 

(21) 

v = We2 /me2. 

To  = H +0$  + p^)/m2c2]1/2, 

0±o  = Piohomc' 

21  Piohom’ 

T = ^v*o- 

A = (to  - -sfte/y0)r, 

A’  = (to  + -sne/y0)r, 

the  functions  and  P4  are  similarly  defined  as  Px  and  P2,  respectively,  but  with  to  replaced 
by  —to,  A replaced  by  —A  , and  A replaced  by  —A.  Finally,  the  double-argument  functions 
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Hs  and  Qs  in  Eqs.  (20)  and  (21)  are  defined  as 

Hs  (a0,  aL ) = Js  ( aL  )ls  (a0,  aL ), 

Qs(a0.  aL)  ==  2HS  (a0,  aL ) + aLJ^(.aL)Is(a0,  aL) 

+ - aLJs  (aL)  (/,_!  (a0,  aL)  — Is+ j (a0,  a^)), 


where 


Is(a0,aL)=  - — f </a7,  (a)a  sin  0 l(a2  - a,2 ) (a22  -a2)]  1/2 


oo 

z 


(at)yj-(a)  cos  s' 


a\  — lflo  azJ» 
a2  — ao  + 

<t>  = sin  [(a2  + al  - a$)H2aaL)]. 


(23) 


(24) 


(25) 


The  integral  series  in  Eq.  (24)  has  been  evaluated  in  Appendix  A with  the  result 


Ma0<  aO  ~Js  (a0)y,'(at). 


(26) 


Substituting  Eq.  (26)  into  Eqs.  (23)  and  (24),  we  obtain 

Hs(aQ,aL)  - (o0  )^j  (flz.)  1 2> 

(27) 

Qs  (fl0>  aL  ) “ 2^5  ) + aLJs  )«C  («/.)Us  (flo  ) (1  + S2/ao  > 

+ \JS  (a0 ) ] 2)  +2  s2Js(a0)Js(a0)Js(aL)laLJs(aL)  - Js(aL)]/a0aL. 

We  assume  that  in  actual  operations,  radiation  at  a particular  cyclotron  frequency  is 
favored,  i.e.  the  parameters  will  be  so  chosen  that  only  one  term  on  the  right  hand  side  of  Eq. 
(20)  or  (21)  dominates.  Hence  we  will  drop  the  summation  signs  in  Eqs.  (20)  and  (21). 
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B.  Physical  interpretation  of  the  beam-wave  coupling  coefficient  Hs  (,knr0,  knrL) 

In  Eq.  (20)  or  (21),  the  first  term  (i.e.  the  term  proportional  to  02o  and  Hs)  originates 
from  the  transverse  motion  of  the  electrons.  This  is  the  source  term  for  the  cyclotron  maser  ra- 
diation. It  can  be  either  positive  (beam  power  gain)  or  negative  (beam  power  loss)  depending 
primarily  on  the  value  of  the  phase  factor  A.  The  second  term  (proportional  to  Qs ) originates 
from  wave  induced  oscillations.  It  is  a positive  term,  hence  leads  to  beam  power  gain.  However, 
unless  the  beam  transverse  velocity  is  too  low  to  be  of  practical  interest,  it  is  always 
' insignificant  when  compared  with  the  first  term.  Thus,  for  practical  purposes,  the  beam  power 

gain  is  proportional  to  Hs  (knr0,  knrL),  which  will  henceforth  be  referred  to  as  the  beam-wave 
coupling  coefficient.  The  physical  significance  of  Hs  can  be  understood  as  follows.  Consider 
the  projection  of  an  electron  orbit  on  the  cross-sectional  plane  of  the  cavity  as  shown  in  Fig.  4, 
* t where  point  Ois  the  axis  of  symmetry,  point  eis  the  position  of  an  arbitrarily  chosen  electron, 

point  C is  the  guiding  center  of  this  electron,  and  the  circle  of  radius  rL  is  its  cyclotron  orbit. 
The  electron  e loses  or  gains  energy  in  the  cavity  fields  through  its  interaction  with  E the 
component  of  the  cavity  electric  field  tangential  to  the  electron  cyclotron  orbit,  i.e. 

E * - Ee  cos  a 
- Eg$J , (k„r)  cos  a, 

(28) 

where  the  z-dependence  of  £ $ has  been  neglected.  We  may  express  E 4 in  terms  of  r0,  rL , and 
<t>  through  the  following  geometrical  relations  (see  Fig.  4): 

r “ (fQ  + r[  — 2 rQrL  cos  0) 1/2 


T 


n 


: 


I! 
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E4  ” E00Jl  t kn  <r0  + rl  ~ 2rOrL  cos  '/21  (rL  ~ r0  cos 

• H + d - 2rorL cos  ~m 

Eeokn~l  'o  l*„  0$  + d ~ 2ro rL  cos  1/21- 


(29) 


Expanding  £ ^ in  terms  of  the  sinusoidal  harmonics  of  <f>  and  noting  that  £ ^ is  an  even  func- 
tion of  0,  we  obtain 


where 


E*  “ L £«rcos 

s -0 


0 7* 

- — J cos  sd> 


(30) 


® V)  ( V ) 1 J d<l>J0  [kn  (r$  + d ~ 2rorL  cos  <*» ) 1/21. 

orL  o 


(31) 


0 


1,  s ■»  0 

2,  s/0. 


Using  tabulated  integral  formulae,44  we  obtain  from  Eq.  (31), 

Eeff SE90  Js  < Vo  < Vi  >• 

SE^H}'2  ( VoVl ) 


(32) 


The  coefficient  £^-is  the  amplitude  of  the  s-th  harmonic  component  of  the  wave  electric 
field  in  the  direction  of  the  electron  velocity.  This  is  the  component  which  provides  the 
effective  electric  field  for  beam-wave  coupling  at  the  s-th  cyclotron  harmonic  frequency.  The 
beam  power  gain  scales  with  the  square  of  Elff>  which  explains  why  it  is  proportional  to 
Hs(knr0,  knrL).  One  notes  that  even  if  the  actual  electric  field  ( EM ) is  large,  the  effective 
electric  field  (£^-)  for  higher  harmonic  interactions  (s  > 1)  may  be  small  because  it  is  propor- 
tional to  H l/2  ( knr0 , knrL).  From  Eqs.  (27)  and  (32),  it  is  easily  seen  that  Eseff  for  s > 1 actually 
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vanished  as  rL  — 0.  This  is  expected  because  non-fundamental  cyclotron  harmonic  interac- 
tions result  from  the  finite  Larmor  radius  effect. 

C.  Normalization 

It  is  convenient  to  introduce  a normalization  scheme  by  which  the  cavity  radius  rw  is 
scaled  out  of  the  results.  This  can  be  achieved  through  the  following  procedures  (normalized 
notations  are  denoted  by  a bar): 

(i)  length  normalized  to  rw  (e.g.  70  - r0/rw)\ 

(ii)  frequency  normalized  to  c/rw  (e.g.  57  — a>rjc)\ 

(iii)  momentum  normalized  to  me  (e.g.  pl0  — p±0/mc)\  and 

(iv)  EM  fields  normalized  to  mc2/erw  (e.g.  — E^erjmc1). 

Other  quantities  such  as  kz  and  r are  to  be  normalized  consistently  with  the  preceeding  pro- 
cedures (e.g.  kz  — kzrw,¥  — r c/rw).  However,  naturally  dimensionless  quantities  such  as 
■y0,  v,  A,  Hs  and  Qs  will  remain  unchanged. 

D.  Optimization  of  guiding  center  position 


In  the  beam  power  gain  function  P (Eq.  (19)]  one  finds  that  the  guiding  center  position 
7q  only  appears  is  the  beam-wave  coupling  coefficient  Hs(.xnr0,  xnrL),  where  the  arguments 
x„70  and  xn?L  are,  respectively,  the  normalized  forms  of  knr0  and  knrL.  Therefore,  to  maxim- 
ize P with  respect  to  70,  one'only  needs  to  find  the  value(s)  of  70  which  maximizes  Hs.  We 
now  examine  the  behavior  of  Hs  ( xnT0 , xnrL)  as  a function  of  its  arguments.  Figures  5a,  b,  c, 
and  d show  the  surface  plots  of  Hs  versus  xn70  and  xnrL  for  the  first  four  cylotron  harmonics. 
We  note  that  the  regime  xnrL  > s is  of  no  interest  because  the  condition  of  frequency  syn- 
chrony (<u  - kzvz0  =»  rfi,/y0)  is  not  satisfiable  there.  From  Fig.  5 and  Eq.  (27),  we  observe 
the  following  properties  of  Hs:  (i)  the  height  of  Hs  dereases  as  s increases;  (ii)  Hs  is  a decreas- 


ing function  of  x„rL  for  s - 1 and  an  increasing  function  of  xnrL  for  s > 1;  and  (iii)  Hs  is  a 
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nonmonotonic  function  of  x„70  with  an  infinite  number  of  peaks.  The  first  peak  is  the  largest 
and  the  subsequent  ones  are  successively  smaller.  Note  that  the  peaks  of  Hs  do  not  in  general 
coincide  with  those  of  Ee  except  for  s — 1.  Given  a fixed  value  of  n,  there  are  values  of  70  for 
which  Hs  vanishes.  Evidently,  it  is  important  to  select  70  such  that  Hs  falls  on  or  near  one  of 
its  peaks  with  respect  to  x„70.  Tables  I,  II,  and  III  list  the  values  of  70  corresponding  to  the 
first  three  peaks  of  Hs  for  various  combinations  of  j and  n.  For  a given  combinatin  of  s and  n, 
the  value  of  70  listed  on  Table  I is  more  preferable  to  those  on  Tables  II  and  III  because  it 
corresponds  to  the  first  peak  of  Hs. 

E.  Total  beam  energy  gain  during  one  transit  time 


The  total  stored  field  energy  ( Wf ) in  the  cavity  is  given  by 

“4>  rl  4 (x„)L/l6. 

(33) 

It  is  convenient  to  introduce  a dimensionless  quantity  F defined  as  the  ratio  of  the  total  beam 
energy  gain  during  one  transit  time  to  the  total  stored  field  energy,  i.e. 


F = Pr/Wf. 


substituting  Eqs.  (19)  and  (33)  into  Eq.  (34),  we  obtain 

F — 2 i/t2  ly0  Jo  (■*«)  "1  _l  (“l  + a2  + “3  + “4). 

where 


(34) 


(35) 


a,  - - Hs  p10  t A ~3  l (w2  - k2  ) (4  sin  2 A/2  - A sin  A) 

+ kz  L A (A  sin  A — 2 sin  2 A/2 ) ] 

+ 2 Qs  A _2  (5T  - k2p2 (,)  sin2  A/2, 

a 2 0 r A "2  A"1  |(«2  — k2  ) [2  (1  + A/A ')  sin  2 A 72  -AsinA'l 


(36a) 
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If  we  decompose  the  cavity  standing  wave  into  a forward  travelling  wave  and  a backward 
travelling  wave,  then  in  Eq.  (8),/£*  and  / L1*  are  the  beam  perturbations  due  to  the  forward 
and  backward  travelling  waves,  respectively.  Similarly  in  Eq.  (35),  a,,  a2  are  due  to  the  in- 
teractions of  /£  * with  the  forward  and  backward  travelling  waves,  respectively,  and  a3  , a4 
are  due  to  the  interactions  of/!iP  with  the  backward  and  forward  travelling  waves,  respective- 
ly. 

In  order  for  the  beam  to  generate,  rather  than  absorb,  electromagnetic  radiation,  it  is  im- 
portant to  match  parameters  so  that  F assumes  negative  values.  The  sign  of  F is  principally 
determined  by  the  phase  factor  A.  Figure  6 shows  a typical  plot  of  Fas  a function  of  A (solid 
curve,  marked  by  F).  The  four  components  of  F represented  by  the  four  terms  in  Eq.  (35)  are 
also  ploted  (dashed  curves  marked  by  1,2,  3,  and  4).  It  is  seen  from  the  plot  of  Fthat  nega- 
tive beam  energy  gain  occurs  when 


— n/7  < A <2 w. 


(37) 
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i 

>■ 


Eq.  (37)  gives  the  resonant  wave  frequency  width  (8w)  and  cyclotron  frequency  width 
(817  e/y0)  for  negative  beam  energy  gain,  i.e. 


8<u  — 2 tt/t, 


(38a) 


and 


8Oe/y0  — 2w/sr. 

(38b) 

One  notes,  however,  that  the  frequency  band  width  of  a cavity  eigenmode  (see  item  F of  this 
section)  is  usually  much  narrower  than  that  given  by  Eq.  (38a). 


Equation  (35)  has  been  evaluated  for  a wide  range  of  parameters,  all  showing  the  same 
condition  [Eq.  (37)  or  (38)]  for  negative  beam  energy  gain.  However,  the  maximum  negative 
value  of  F (point  X in  Fig.  6)  shifts  toward  smaller  A when  L increases.  Also,  F increases  shar- 
ply when  r and  P |0  increase,  as  expected  from  the  analytical  form  of  F.  Here  we  note  an  im- 
portant difference  between  the  cyclotron  maser  interactions  in  waveguides  and  cavities.  In  a 
waveguide,  only  one  type  of  interaction  is  present,  namely,  the  beam  and  forward  travelling 
wave  interaction.  It  has  been  shown  1 ~3  that,  for  such  a case,  negative  beam  energy  gain  oc- 
curs only  when  w ~ kzPzO  ~ s^f/yo  ^ “A/t)  >0.  This  interaction  corresponds  to  the  a j 
term  of  Eq.  (35).  As  shown  in  Fig.  6,  the  a,  term  indeed  assumes  negative  values  only  when 
A > 0 (consider  only  the  region  — 2ir  < A < 2n  ).  In  a cavity,  the  saturation  is  more  com- 
plicated because  the  beam  interacts  with  a standing  wave  (two  oppositely  directed  travelling 
waves).  Consequently,  there  are  four  types  of  interactions  leading  to  the  four  terms  in  Eq. 
(35),  as  just  pointed  out.  Furthermore,  all  four  terms  are  important  (see  Fig.  6)  and  when  they 
add  up,  the  net  effect  is  such  that  negative  beam  energy  gain  occurs  not  only  for  positive  A but 
for  slightly  negative  A also  (see  Fig.  6).  The  differences  in  the  interaction  processes  may  also 
have  important  implications  in  nonlinear  considerations.  For  example,  in  order  to  enhance  the 
beam-to-wave  energy  transfer  efficiency  in  a waveguide,  one  only  has  to  raise  the  saturation 
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level  of  the  beam  and  forward  travelling  wave  interaction  through,  for  example,  magnetic  field 
adjustment 13,14  . In  a cavity,  however,  raising  the  saturation  level  of  one  type  of  interaction 
does  not  necessarily  lead  to  a higher  overall  efficiency  since  there  are  three  other  types  of  in- 
teractions whose  contributions  might  be  adversely  affected. 


In  view  of  the  lengthy  algebra  involved  in  the  present  analysis,  an  independent  check  of 
the  results  appears  warranted.  We  have  developed  a numerical  code  which  calculates  a single 
electron  orbit  (hence  its  energy  change)  in  the  cavity  fields  of  Eqs.  (2)  to  (4).  The  field  ampli- 
tude is  kept  sufficiently  low  so  that  orbit  perturbations  remain  linear.  Taking  an  ensemble 
average  over  the  distribution  represented  by  Eqs.  (17)  and  (18),  we  then  obtain  the 
corresponding  values  of  Ffrom  the  numerical  code,  which  are  shown  by  dots  in  Fig.  6.  The 
agreement  between  the  numerical  and  analytical  approaches  has  been  excellent  for  the  case 
shown  and  also  for  a large  number  of  other  cases  we  have  checked. 

F.  Threshold  beam  power 

The  present  model  assumes  an  idealized  cavity,  namely,  a cavity  with  discrete  eigenfre- 
quencies.  In  practice,  a cavity  always  has  a finite  Q [defined  in  Eq.  (39)]  due  to  loading  or  wall 
resistivity,  which  leads  to  narrow  bands  of  eigenfrequencies  with  width  ~ 2 u>/Q.  From  Eq. 
(38a),  we  find  that  the  frequency  width  for  negative  beam  energy  gain  is  approximately  2ir/r. 
Thus,  a cavity  of  finite  Qcan  be  effectively  treated  as  an  idealized  cavity  if 

7t/t  :»  iti/Q. 

This  is  a condition  easily  satisfiable  even  for  unusually  low  values  of  Q or  large  values  of  t. 


Because  of  finite  Q,  energy  is  drained  from  the  cavity  at  the  rate 

while  the  beam  pumps  energy  into  the  cavity  at  the  rate 


(39) 
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From  the  condition  of  power  balance,  Pjn  ^ 
oscillations, 


- FWj! T. 


(40) 


Poun  we  obtain  a threshold  condition  for  cavity 


The  beam  power  is  given  by 


— FQ^  <DT 


(41) 


pb  “ H(yo  -Dmc2*#. 

Inserting  Eq.  (35)  into  Eq.  (41)  and  making  use  of  Eq.  (42),  we  may  rewrite  Eq.  (41)  as 


(42) 


P„  > P 


Th 

b 


where  P™,  the  threshold  beam  power  for  cavity  oscillations,  is  defined  as 

tv,  _ o (yp  -1)-/qUb)  m2cs  ^ 

b QL(an  + a2  + «3  +<*4)  2 x 1010e2 


(44) 


In  Eq.  (44),  c,  e,  and  m are  to  be  expressed  in  Gaussian  units,  while  all  other  quantities  are  ei- 
ther normalized  or  dimensionless. 


Equation  (44)  is  a function  of  nine  free  parameters,  i.e.  s,  n,  /,  70,  L,  y0,  Q,  ft,,  and  the 
ratio  of  electron  transverse  velocity  to  parallel  velocity  (v^/v^).  In  the  numerical  data  to  be 
presented,  we  let  the  radial  eigenmode  number  n equal  the  cyclotron  harmonic  number  s.  Such 
a combination  of  n and  s allows  us  to  choose  an  optimum  r0  from  Table  I (see  the  numbers 
marked  by  *)  which  is  approximately  halfway  between  the  cavity  wall  and  the  axis.  The  axial 
eigenmode  number  / will  be  fixed  at  / — 1.  With  s,  n,  and  / known,  the  cavity  eigenfrequency 
to  can  be  calculated.  The  cyclotron  frequency  ft  , will  be  so  specified  that  the  beam  energy  loss 
is  a maximum  with  respect  to  A (see  point  X on  Fig.  6).  The  velocity  ratio  will  be  fixed  at 
vjVvzO  ” 1-5,  which  is  consistent  with  most  experimental  conditions.  Finally,  Q and  Pjh  can 
be  combined  into  a single  parameter.  The  remaining  parameters  s,  L , and  y0  will  be  varied. 
Note  that  the  cavity  wall  radius  rw  has  been  scaled  out  of  the  final  results.  Figures  7a,  b,  and  c 
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plot  QPbh  as  a function  of  the  electron  kinetic  energy  Wb,  where  Wb  — (y0  - 1 )mc2,  for 
the  first  four  cyclotron  harmonics  and  for  L — 5,  10,  and  20.  It  is  seen  that,  for  a fixed  Q , the 
first  cyclotron  harmonic  can  be  excited  with  the  lowest  threshold  beam  power  Pbh,  and  the 
higher  the  beam  energy,  the  higher  Pbh  becomes.  For  higher  cyclotron  harmonics,  on  the 
other  hand,  Pbh  is  considerably  higher,  but  there  is  an  optimum  beam  energy  for  which  QPbh 
is  lowest.  Although  P[h  can  be  lowered  by  increasing  either  Q or  L,  for  reasonable  values  of  Q 
and  I,  it  still  requires  a relatively  high  electron  energy  and  beam  power  to  excite  cyclotron  har- 
monics above  2.  As  an  example,  we  apply  the  present  results  to  the  experiments  of  Kisel  et. 
al.35  Using  their  parameters  (Q=  4000,  L = 7,  and  Wb  = 20  keV),  we  find  from  Eq.  (44)  [or 
approximately  by  interpolating  from  Figs.  7(a)  and  (b)l  that  Pbh  for  the  first,  second,  third  and 
fourth  cyclotron  harmonics  are,  respectively,  1.5  kW,  22.1  kW,  297.5  kW,  and  3800  kW.  Thus, 
their  beam  power  (<  70  kW)  was  sufficient  to  excite  the  first  and  second  cylotron  harmonics 
(as  reported  in  their  paper)  but  far  below  the  necessary  power  to  excite  the  third  and  fourth 
cyclotron  harmonics.  It  is  also  found  from  Eq.  (44)  that  their  beam  energy  (—20  keV)  is  very 
close  to  the  optimum  energy  (i.e.  the  energy  that  requires  the  lowest  threshold  beam  power) 
for  exciting  the  second  harmonic  but  is  a factor  of  3 and  8 below  the  optimum  energies  for  the 
third  and  fourth  cyclotron  harmonics,  respectively. 
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V.  DISCUSSION 


We  first  consider  the  validity  of  the  linear  approximation.  In  solving  the  Vlasov  equation, 

} 

we  have  made  the  assumption  (in  Gaussian  units) 

^ So 

(45) 

Making  use  of  Eqs.  (33)  and  (39),  this  assumption  can  be  rewritten  as  a condition  on  the  out- 
put power  Pout 

QPou,  « TT  fib2  rl  Jo  (x„)mL. 

16  (46) 

where  Pow  is  the  sum  of  the  output  wave  power  and  the  power  dissipated  in  walls.  This  condi- 
tion is  generally  well  satisfied  at  the  threshold  beam  power  level.  As  a numerical  example,  we 
use  the  parameters  of  Fig.  6 and  assume  rw  — 1 cm,  then  Eq.  (46)  gives 

QP0U,  « 3x10 7 

(47) 

For  the  same  parameters,  QP Jh  is  only  5x10 4 kW  (see  Fig.  7a). 

Using  Eq.  (37),  we  may  derive  a validity  condition  for  the  cold  beam  assumption.  Equ  - 
tion  (37)  shows  that  the  range  of  A for  negative  beam  energy  gain  is  approximatly  2.5ir,  where 


A has  been  defined  as 


Thus,  if 


A - (<*>  ~ Mrf  ~ sn/yoJF. 


(kz  8/3 z 4-  sfl,,  8y/y^)  r <3C  2.57T, 


where  hpz  and  by  are,  respectively,  the  axial  velocity  spread  and  energy  spread,  the  cold  beam 
assumption  is  essentially  valid.  Since  kz  ” lir/L  and  r “ we  may  rewrite  condition  (49) 
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(/tt8/32  +s(leL8y/y q)/^^  ^ 2.5ir. 


Thus,  the  condition  becomes  more  restrictive  as  the  axial  eigenmode  number  (/),  the  cavity 
length  (£),  and  the  wave  frequency  (—  sCle/y0 ) increase  or  as  the  beam  axial  velocity  (P^) 
decreases. 


Turning  to  nonlinear  considerations,  we  can  obtain  a crude  estimate  of  the  limiting  cavity 
length  beyond  which  the  energy  conversion  efficiency  begins  to  drop.  Equation  (37)  gives  the 
approximate  range  of  A for  which  the  beam  will  lose  energy  in  traversing  the  cavity,  where  in 
the  linear  analysis  A has  been  defined  in  Eq.  (48).  Since  Eq.  (37)  is  a phase  synchrony  condi- 
tion, we  expect  it  to  hold  true  in  the  nonlinear  stages  of  the  interaction  as  well.  In  which  case 
we  replace  the  initial  energy  -y0  'n  Ed-  (48)  with  the  instantaneous  energy  y and  note  that  p^ 
remains  relatively  constant  because  of  the  absence  of  any  axial  electric  field.  Assuming 
< y f>  is  the  average  electron  energy  after  leaving  the  cavity,  then  the  average  change  of  A 
due  to  the  change  of  y is  approximately  ( < y f>  -1  - -y0_1 ).  Thus  condition  (37)  can 
be  written 

sfler(y0  - < yf>  )/y0  < y r>  < 2.5ir 

(51) 

Letting  17  = (yQ  — <yy>  )/(-y0  — 1 ) be  the  energy  conversion  efficiency,  and  writting 
L = P^T,  we  obtain  from  Eq.  (51) 

Y<  y _ 37rffzQ>,ojyo  -vhp  ~1)1 
siler)  ( y0  — 1 ) 


Equation  (52)  gives  the  limiting  cavity  length  Lc  for  achieving  maximum  efficiency.  Note  that 
Lc  decreases  as  the  efficiency  and  wave  frequency  increase  or  as  the  beam  axial  velocity  de- 
creases. Thus,  while  it  is  possible  to  lower  the  threshold  beam  power  by  increasing  the  cavity 
length,  this  may  also  lower  the  efficiency  if  condition  (52)  is  violated. 
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To  carry  out  the  integration  in  Eq.  (A.4),  we  replace  the  variable  of  integration  a with  x,  where 
x is  defined  through  the  equation 
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Table  I — Values  of  r0  for  which 
H,  (xnr0,  xnrL)  falls  on  its  first 
peak  with  respect  to  the  argument 
x„Tq.  The  numbers  marked  by  * 
have  been  used  to  obtain  the  data 
presented  in  Pig.  7. 


\ s 

l\ 

1 

2 

3 

4 

1 

0.48* 

0.80 

2 

0.26 

0.43* 

0.60 

0.76 

3 

0.18 

0.30 

0.41* 

0.52 

4 

0.14 

0.23 

0.32 

0.40* 

5 

0.11 

0.19 

0.26 

0.32 

6 

0.09 

0.16 

0.21 

0.27 

7 

0.08 

0.13 

0.18 

0.23 

8 

0.07 

0.12 

0.16 

0.21 

Table  II—  Values  of  rQ  for  which 

Hi  (xnT 0*  xnW  fal18  on  its  se‘ 
cond  peak  with  respect  to  the 

argument  xnT0. 


1 

2 

3 

4 

0.76 

0.87 

0.52 

0.60 

0.79 

0.91 

0.40 

0.46 

0.60 

0.70 

0.32 

0.37 

0.49 

0.56 

0.27 

0.34 

0.41 

0.47 

0.23 

0.29 

0.35 

0.41 

0.21 

0.26 

0.31 

0.36 

Table  III  — Values  of  r0  for  which 
Ht  (xnT0,  xnrL)  falls  on  its  third 
peak  with  respect  to  the  argu- 
ment *nF0. 


X 

1 

2 

3 

4 

1 

2 

3 

0.84 

0.98 

4 

0.64 

0.75 

0.85 

0.95 

5 

0.52 

0.61 

0.69 

0.77 

6 

0.44 

0.51 

0.58 

0.65 

7 

0.38 

0.44 

0.50 

0.56 

8 

0.33 

0.38 
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Fig.  1 - (a)  End  view  of  the  election  cyclotron  mater  configuration. 
Guiding  centers  of  all  electrons  are  uniformly  distributed  on  the  circle 
of  constant  radius  fq.  The  applied  magnetic  field  (not  shown)  points 
toward  the  reader,  (b)  Side  view  of  the  same  configuration. 
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Fig.  6 - F (tolid  wave)  versus  A for  * ■ n « / * 1,  Pq  * 0.48,  T « S,  7q  * 1.1,  and 
vjjyvz0  * 1.5.  The  four  components  of  F - the o “3>  »n<*  °4  terms  in  Eq.  (35)-are 
also  plotted  in  dashed  curves  (marked  by  1,  2,  3,  and  4,  respectively).  Solid  dots  are  the 
values  of  F calculated  from  an  independent  [numerical  code.  Point  marked  by  X is  the 
maximum  negative  value  of  F, 
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Fig.  7 - Threshold  beam  power  ( P ^ ) times  Q versus  beam 
kinetic  energy  (H^for  / ■ 1,  n » *,  vjj/vz0  = 1J,  and  (a)  T 
= 5,  (b)  L = 10,  (c)  L • 20.  The  beam  guiding  center  position 
~q  is  so  chosen  that  Hs(x^q,  x^)  falls  on  its  first  peak  with 
respect  to  (see  Fig.  5).  The  magnetic  field  Bq  has  been 
specified  to  maximize  the  beam  energy  low  (see  point  X 
on  Fig.  6). 


